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Verfahien zum Segmentieren einer dreidiniensionalen Struktur aus einem dieidimensionalen 
Datensatz unter Beriicksichtigung von Benutzeatorrekturca 

Die BjBndung betrifft ein Verfehrm zum Segmentieren einer dieidimensionalen Struktur aus 
5 einem dieidimensionalen, insbesondere medizinischen Datensatz, unter Beriicksichtigung von 
Benutzericorrekture^ Das Veifahren wird mit BBlfe eines dieidimensionalen deformierbaien 
Modells, dessen Oberflache durch ein Netz mit Netzpunkten und diese verbindende Maschea 
gebildet ist, durchgefuhrt Die Erfindung bezieht sich auBerdem auf eine Bildveraibeitungsein- 
richtung zur Durchfuhrung dieses Verfahrens sowie auf ein Computerprogramm zur deren 
10 Steuerung. 

En Verfstfiien der eingangs genannten Art, bei dem aUerdings a priori keine Benutzerkorrek- 
tuien voigesehen sind, ist aus einem Aufsatz von Mclnerney et aL , JDefoimable models in 
medical image analysis: A survey** in der Zeitschrift Medical Image Analysis, 1(2):91- 108, 

15 1996 bekannt Danach wird ein deformierbaies Modell als eine elastische Oberflache dar- 
gestellt, die unter dem Finfluss einer internen und einer external sogenannten Enetgie ihre 
Form und Position vecandem kann. Dabei hat die interne Enagie die Aufgabe, die Form des 
Modells, das aufgrund medizinischen Vorwissens iiber die zu segmentierende Struktur erzeugt 
worden sein kann, moglichst gut zu erhalten, wahiend die exteme Eheigie ein Bewegen der 

20 Mo delloberflache in Richtung auf die Struldnrrander bewirkt Die exteme Energie wird aus 
dem dieidimensionalen Datensatz abgeleitet Die Maschen haben zJ3. die Form eines Drei- 
ecks, dessen Eckpunkte durch drei benachbarte Netzpunkte an der Oberflache des Modells 
defirriert sind. 



25 Bei dem bekannten Verfahien weiden in dem dieidimensionalen Datensatz entlang von 
FMchennormalen der Dreiecksflachen Kandidatenpunkte gesucht, dh. Punkte, die wahr- 
scheinlich auf der Oberflache der zu segmentierenden Struktur angeordnet sind, wobei die 
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Oberfiache in der Regel dutch einen ausgepragten Gradienten der Datenwerte gekennzeichnet 
ist Nachdem fur die Mascben des die Oberfiache des ModeUs definiesrenden Netzes Kandi- 
datenpunkte ermittelt worden sind, erfolgt eine Neuberechnung der Netzpunkte des ModeUs 
auf Basis der gefundenen Kandidatenpunkte. Dabei werden die neuen Netzpunkte so berech- 
5 net, dass fur diese die gewichlete Summe aus interner und extemer Eaergie ein Minimum ein- 
nirnmt Danach wird das Verfahren unter Zugrundelegung der neu berechneten Netzpunkte 
mehrfach wiederholt, wobei sich beijeder Wiederholung das Model der Form der zu seg- 
mentierenden Struktur annahert Dieses iterative Verfahren ist beendet, sobald ein Abbruch- 
kriterium, z3. eine bestimrnte Anzahl an Wiederholungen, erfuUt ist 

10 

Nachteilig ist bei diesem und auch bei anderen automatischen Segmentierungsverfahren, dass 
Segmentierungsfehler entstehen, z.B. durch Artefakte in dem dreidimensionalen Datensatz. 
Diese Segmentierungsfehler treten haufig auf, wenn der Abstand einer zu segmentierenden 
Struktur zu einer anderen Struktur klein ist. Zur Vermeddung von Segrnentierungsfehlem sind 
15 verschiedene Verfahren bekannt, von denen im folgenden einige erlautert werden. 

Ein einfaches Verfahren zur Vermeidung bzw. zur Konektur von Segrnentierungsfehlem ist die 
Wiederholung des oben beschriebenen Anpassungsverfahrens mit unterschiedlichen Anfangs- 
modeUen und unterschiedlichen Anfangspositionen des aeformierbaren ModeUs relativ zu der 
20 zu segmentierenden Struktur. Dieses Verfahren erfordert ein Vorhersehen des Ergebnisses 
des Anpassungsverfahrens und somit viel Erfahrung des Anwenders. Dieses Vorhersehen ist 
nahezu unmoglich fur komplexe zu segmentierende Strukturen. Aufierdem ist dieses Verfahren 
sehr zeitintensfv. 

25 Ein gattungsgemHBes Verfahren gibt dem Anwender dagegen die MOgtichkeit, wahrend des 
Anpassungsverfahrens in den AnpassungsprozeB einzugreifen, indem das ModeU bspw. durch 
Bewegen von einem oder mebreren Netzpunkten manuell verformt wird Nachteilig ist bei 
diesem Verfahren, dass die Visuahaerung von Benutzeikorrekturen nicht in iibersichtHcher 
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Weise fur die gesamte zu segmentierende Struktur moglich ist Deshalb muss das deformier- 
bare Modell schichtweise kontrolliert und gegebenenfalls konigiert werden. Eine genaue Kor- 
rektur ist daher sehr zeitintensiv. 

5 Bin weiteres bekanntes, gattungsgemaBes Verfahren besteht daxin, dass, nachdem eine auto- 
matische Segmentierung stattgefunden hat, das angepasste Modell, z3. durch Verschieben 
eines Netzpunktes, manuell verforait wind Mit diesem verformten Modell wild dann die auto- 
matische Segmentiexung ein zweites Mai durchgefuhit Problematisch ist bei diesem Verfah- 
ren, dass der Verfahiensschritt, der die interne und die exteme Energie minimiert, die manuell 
10 verschobenen Netzpunkte zuriick an ihie urspriinglichen Positionen bewegt, da diese der minv 
malen Energie des defonrrierbaren Modells entsprechen. 

Es ist Aufgabe der vorliegenden Erfindung, ein Verfahren anzugeben, bei dem Segmentfe- 
rungsfehler mit einem geringen Zeit- und Rechenaufwand vermieden bzw. konigiert werden 
15 konnen. Diese Aufgabe wild erfindungsgemaB gelost durch ein Verfahien zum Segmentieren 
einer dreidimensionalen Struktur aus einem dreidimensionalen, insbesondere medizinischen 
Datensatz, unter Berucksichtigung von Benutzerkorrekturcn mit den Schritten: 

a) Bereitstellen eines dreidimensionalen, deformierbaren Modells, dessen Oberflache durch 
20 ein Netz mit Maschen gebildet ist, die Netzpunkte an der Oberflache des Modells 

miteinander verbinden, 

b) Anordnen des Modells an einer Stelle im drcidiiiiensionalen Datensatz, an der sich die zu 

segmentierende Struktur befindet, 

25 

c) manuelles Verschieben von Netzpunkten, 

d) Neuberechnen der Netzpunkte des Modells unter gewichteter Berucksichtigung von 

manuell verschobenen Netzpunkten. 
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Im Gegeiisatz zu bekannten Verfahren werden bei der Erfindung die manuell verschobenen 
Netzpunkte bei der Neuberechnung des deformierbaren Modells gewichtet berucksichtigt 
Dies ermSglicht eine untexschiedlich starke Beriicksichtigung jedes manuell verschobenen 
5 Netzpunktes bei der Deformation, was zu einer Benutzerkonektur mit geringem Zeitaufwand 
fuhrt 

Ein bevorzugtes Verfahren, bei dem verschobene Netzpunkte, deren Abstand zu einer Grenz- 
flache der zu segmentierenden Struktur klein ist, starker berucksichtigt werden als verschobe- 
10 ne Netzpunkte mit einem grQBeren Abstand zu dieser Grenzflache, beschreibt Anspruch 2. 

Arepruch3 definiert embevorzugfcesVerfaih^ daszu 
guten Segrnentiennigpergebrtt^^ 

15 Eine Bild\osrbdtun^aririchrung zur Durchfuhrung des edSndun^^mafim Verfaihrens ist 
in An^ruch 4 beschrieben. Atispruch 5 definiert ein Connputerprograrnm zur Steuerung der 
Bildveraibdtun^eirinchh^ nach Aispruch 4. 

Die Erfindung wird im folgenden anhand der Zeichnungen naher erlautert Es zeigen: 

20 

Fig. 1 eine zur Ausfuhrung des erfindungsgemaBen Verfahrens geeignete 
BUdverarbtitragseiimchtung in schematischer Darstellung, 
Fig. 2 ein Ablauf diagramm des erfindungsgemaBen Verfahrens, 
Fig. 3 ein deformierbares Model! eines Riickenwirbels, 
25 Fig. 4 einen Teil einer Wirbelsaule mit einem auf einem Ruckenwirbel angeordneten 
deformierbaren Model], 

Fig. 5 das auf dem Ruckenwirbel angeordnete deformierbare ModeH nach Ausfuhrung eines 
bekannten Anpassungsverfahrens, 
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Fig. 6 das auf dem Riickenwirbel angeordnete defonnierbare Modell nach einer manuellen 
Korrektur und 

Fig. 7 das auf dem Riickenwirbel angeordnete defonnierbare Modell nach Ausfuhrung des 
erfindungsgemaBen Verfahiens. 

5 

Die in Fig. 1 daigestellte Bildve^bdtxmgseiiirichtung umfasst einen Bildveraibeitungs- und 
Steuerprozessor 1 mit einem Speicher 2, in dem ein deformierbares Modell einer zu seg- 
roentierenden Struktur und ein dreidimensionaler, insbesondere medizinischer Datensatz ge- 
speichert sein konnen. Der Bildverarbeitungs- und Steuerprozessor 1 karm z.B. mittels eines 

10 Glasfaserkabels 3 mit einer nicht naher dargestellten bildgebenden Einrichtung, wie einem 
Magnetresonanz- oder einem Computertomographiegerat, gekoppelt sein. Die von dem 
Bildveraibeitungs- und Steuerprozessor 1 segmentierte Struktur kann auf einem Monitor 4 
wiedergpgeben werden. Der Benutzer kann uber eine Tastatur 5 oder durch andere, in Fig. 1 
nicht naher daigestellte Eingabeeinheiten auf den Bildverarbeitungs- und Steuerprozessor 1 

15 zugreifen. 

Fig. 2 zeigt den Ablauf eines Segmentierungsverfahrens, das mit der Bildveraibeitungs ein- 
richtung nach Figur 1 durchgefuhrt werden kann. 

20 Nach der Ihitialisierung im Schritt 100 wird im Schritt 101 das defonnierbare Modell M einer 
zu segmentierenden Struktur sowie ein drcidimensionaler Datensatz geladen, der die zu seg- 
mentierende Struktur enthalt 

Fig. 3 zeigt ein deformierbares Modell M, das in diesem Ausfuhrungsbeispiel das Modell eines 
25 Rixckenwirbels (vertebra) ist Auch andere Struktuien kOnnen mit BBlfe eines defo mrierbaren 
Modells segmentiert weiden, z.B. der am oberen Bide des Oberschenkelknochens befind- 
HcheFemuikopf oder ein Organ, bspw. das Herz. Die Oberflache des deformierbaien Mo- 
dells M wird durch ein Netz gebildet, dessen dreieckffcrmige Maschenje drei benachbarte 
Netzpunkte auf der Oberflache des Modells miteinander verbinden. Anstelle eines in dieser 
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Weise triangulierten Modells konnte die OberMche auch durch eine beliebige Polygonstruktur 
definiert sein, bei der die Maschen des Netzes nicht dreieckformig sind, sondern eine von drei 
verschiedene Anzahl an Netzpunkten unter Bildung eines Polygons miteinander verbinden. 

5 Ein deformierbares Modell kann aus mehreren Trainingsdatensatzen erzeugt werden, in denen 
jeweils das gleiche Objekt dutch einen Fachmann so manuell segmentiert worden ist, dass 
jedes segmentierte Objekt die gleiche Topologie aufweist Eine Mittelung, im einfachstea FaU 
eine arithmetische Mittelung, tiber alle segmentierten Objekte fuhrt zu einem Objekt, dass als 
deformierbares Modell bezeichnet wird. 

10 

Das erfindungsgemaBe Verfahren ist auch anwendbar, wenn das deformierbare Modell in 
anderen Ausfuhrungsformen auf andeie Art und Weise erzeugt worden ist, zB. durch 
Tnangulation eines Modells eines Ruckenwirbels. Ein haufig angewandtes Triangulations- 
verfabren ist aus dem Artikel .Marching cubes: A high resolution 3d surface construction 
15 algorithm.", WJE. Lorensen und H£. dine, Computer Graphics, 21(3):163-169, 1987 
bekannt. 

hn Schritt 102 wird das deformierbare Modell M relativ zur zu segmentierenden Struktur so 
ausgerichtet, dass die Oberflache des Modells mogUchst gut mit den Konturen der zu seg- 
20 mentierenden Struktur UbereinstimmL 

Tefle des deformierbaren Modells konnen im Schritt 103 manuell verandert werden. Dabei 
karm ein Anwender Dreiecksflachen oder einzelne Netzpunkte des Modells verschieben. hn 
Schritt 104 wird furjede Masche bzw. in diesem Austubrungsbeispiel fur jedes Dreieck der 
25 Oberflache des deformierbaren Modells jeweils ein Kandidatenpunkt x, ermittelL Der Kan- 
didatenpunkt x, liegt auf derHMchennormalen ik eines Dreiecks, wobeidie Flachennormale 
n ein Einheitsvektor ist und von einem Mttelpunkt Xi der Dreiecksflache ausgeht Die 



-7- 



FHDE020307EP-P 



Positionen der Mittelpunkte bezeichnen die Positionen der entsprechenden Dreiecksflachen 
bevor sie im Schiitt 103 u JU. manueJl verschoben worden sind Die Position des Kandidaten- 
puriktes auf der Flachennormalen wird nun so gewahlt, dass die Funktion 

fU) = Dj 2 8 2 - F(x, + j8n,) (1) 

5 

ein Minimum aufweist Hierbei ist D ein frei wahlbarer, aber konstanter Wert, der angjbt, in 
welchem AusmaB der Abstand zur Drciecksflache bei der Bmitflung eines Kandidaten- 
punktes beriicksichtigt wind. Des weiteren ist d ein dem Abstand zweier auf der FBchennor- 
malen aufeinander folgender Punkte entsprechender Skalar und j ein Index, der Punkte auf 
10 der Flachennormalen bezeichnet Die sogenannte KantendeteMonsfunkfion F(x) ist ein MaB 
fiir die Wahrscheinlichkeit, dass der Punkt x auf der Oberflache der zu segmentierenden 
Struktur liegL Diese Funktion kann bspw. nach der Beziehung 

F(x)=±n r g(x) < 2 > 

15 berechnet werdea In dieser Gleichung ist g(x) der Gradient der Datenwerte, dies sind in 

diesem Ausfuhrungsbedspiel Helligkeitswerte, des dreidimensionalen Datensatzes an der Stelle 
x. Da die Funktion F(x) dem Skalarprodukt der Vektoren n und g(x) entspricht, wird die 
Funktion besonders groB, wenn der Gradient in Richtung der FlMiennormalen verlauft DJl, 
wenn angenommen wird, dass sich an der Grenzflache der zu segmentierenden Struktur die 

20 Helligkeitswerte des dieidimensionalen Datensatzes stark andern, die Funktion F(x) 
besonders groB ist, wenn die Hachennormale senkrecht auf dieser Grenzflache steht In 
diesem Ausffihrungsbeispiel wird angenommen, dass im Bereich der Gr enzfl a che die zu 
segmentierende Struktur Datenwerte aufweist, die grOfier sind als Datenwerte auBerhalb 
dieser Struktur, Es wird also angenommen, dass der Gradient g(x) in die zu segmentierende 

25 Struktur weist Die HSchennormale n weist dagegen nach auBen. Das Vorzeichen der 
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Funktion F(x) in Gleichung (2) ist in diesem Ausfuhrungsbeispiel nun so gewahlt, dass im Falle 
einer wie hier definierten Grenzflache die Funktion F(x) positiv wird In anderen Ausfuhrungs- 
beispieien kann das Vorzeichen auch so gewahlt werden, dass bei den oben genannten 
GroBenverhaltaissen der Datenwerte im Bereich der Gienzflache die Funktion F(x) negativ ist 

5 

Der Kandidatenpunkt x, berechnet sich nach folgendear Gleichung: * 

Ji*.- (3) 

Hierbei bezeichnetj* denjenigen ganzzahligen Wert von j, fur den die Funktion der Gleichung 
(1) ein Minimum aufweist 

10 

Nachdem Dreiecksflachen oder Netzpunkte manuell verschoben worden sind und fur jede 
Dreiecksflache der Oberflache des deformierbaien Modells jeweils ein Kandidatenpunkt 
ermittelt worden ist, erfolgt im Schritt 105 eine Anpassung der Oberflache des deformierbaien 
Modells an die gefundenen Kandidatenpunkt© und die manuell verschobenen Dreiecksflachen 
15 oder Netzpunkte. Dabei wird das Modell so deformiert, dass die sogenannte Energiefunktion 

(4) 

+aE int + P£ Bser 

ein Minimum ednnimrnt EBerbei ist E* t die exteme Energie, die die Dreiecksflachen in Richtung 
der im Schritt 104 ermittelten Kandidatenpunkten bewegt, wahrend E nt die interne Biergje ist, 
die einer Deformation des Netzes entgegenwirkt Die Energie ist eine Energie, die eine 
20 Deformation des Modells in Richtung der vom Anwender verschobenen Netzpunkte oder 
Dreiecksflachen bewirkt Die Gewichtungsfaktoren a und 6 sind typischerweise so gewahlt, 
dass die drei Energien im wesentHchen in gleichen Teilen zur Gesamtenergie beitragen. 



25 
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Als exteme Energie kann jede Qrofie herangezogen werden, die sich bei einer Annaherung 
des deformierbaren Modells an die dutch die Kandidatenpunkte definierte Struktur 
veikleinert In diesem Ausfiihrungsbeispiel wild die exteme Energie ^ durch folgende 
Den: 

A 

5 

wobei tiber samtliche Dreiecke summiert wiixL Der Gewichtungsfaktor w t berechnet sich nach 
folgender Gleichung: 

w t = max {),F(x f ) - Dj*8 2 } (6) 

GemaB Gleichung (6) ist der Gewichtungsfaktor fur Kandidatenpunkte x, , die mit einer 
10 groBen Wahrscheinlichkeit auf der Gienzflache der zu segmentieienden Struktur liegen, die 
also einen hohen F(x.) - Wert aufweisen, und deren Abstand zu den Dreiecksflachen gering 
ist, groBer als fur andere Kandidatenpunkte. 

Die interne Energie Ei nt kann jede Grofie sein, die umso kleiner wird, je weniger das Modell 
15 deformiext wird. In diesem Ausfiihrungsbeispiel berechnet sich intone Energie nach der 
Beziehung 

y (7) 
^int = £ £ - ) - sRA,,) 2 . 

Hierbei ist V die Anzahl der Dreiecksflachen des Modells, s tin Skaheningsfokfcor, R eine 
Rotadonsmatrix und A u der Abstand zwischen zwei Mittelpunkten zweiea: Dreiecksflachen 

20 des ursprunglichen, undefcamieiten Modells. Die Menge N(0 beinhaltet die Mittelpunkte 
derjenigen Dreiecksfifichen, die zur Dreiecksflache mit dem Mittelpunkt x { benachbart sind 
Bd der ea^troaHg^Minimieamg derEneigie gemMB Gleichung (4) ist s = l, und R ist die 

Einhe i tematrh r , 
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Die Energie E usa kannjeder Ausdruck sein, der kleiner wird, wenn das deformierbare Modell 
den verschobenen Dreiecksflachen bzw. Netzpunkten angenShert wird, und der die Beitoage 
der einzelnen verschobenen Dreiecksflachen bzw. Netzpunkte wichteL Die Energie be- 
5 rechnet sich in dieser Ausfuhrungsform nach folgender Gleichung: 

£ nse r=Ew t (x t -x < ) 2 , (8) 

A 

wobei tiber alle Dreiecke, die im Schritt 103 manuell verandert worden sind, summiert wird 
und % die Mittelpunkte dieser Dreiecke sind. Die Multiplication mit den Gewichtungsfaktoren 
10 w, = max{ 0, F(x ( ) } bewirkt, dass diejenigen Dreiecksflachen, die von dem Anwender auf 
eine Grenzflache der zu segmentierenden Struktur verschoben worden sind, starker berttck- 
sichtigt werden als andere Dreiecksflachen. 

In anderen Ausfuhrungsformen konnte durch die Gewichtungsfaktoren w t zusatzlich der 
15 Abstand der Mittelpunkte der manuell verschobenen Dreiecksflachen zu den Kandidaten- 
punkten beriicksichtigt werden. Dies konnte bspw. durch die Gewichtungsfaktoren 

(9) 

iv, = maxfa F(x f ) - Dp., - x ,]} 

eneicht werden. Dann wurden manuell verschobene Dreiecksflachen, die in der Nahe der 
Oberflache der zu segmentierenden Struktur platziert worden sind, starker beriicksichtigt 
20 werden als verschobene Dreiecke, die weiter entfernt von der zu segmentierenden Struktur 
sind. 

« 

Nachdem durch Minirnierung der Energie gemSB Gleichung (4) dieMtterpunkteXi der 
Dreiecksflachen des deformierten Modells ermittelt worden sind, werden ein neuer 
25 Skaherungsfaktor s und eine neue Rotationsmatrix R bestirnmL Dabei werden unter 
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Beriicksichtigung der aktuellen Mittelpunkte Xi der Skalierungsfaktor und die Rotationsmatrix 
so gewahlt, dass die interne Energie minimal wird In anderen Ausfuhrungsfonnen konnten die 
Mittelpunkte x i5 der Skalierungsfaktor s und die Rotationsmatrix R gleichzeitig bestimmt wer- 
den, indem bei der Bestimmung des Minimums der Gesamteneargie E audi der Skalierungs- 
5 faktor und die Rotationsmatrix variabel sincL 

Nachdem im Schritt 105 das deformierbare Modell an die zu segmentierende Struktur ange- 
passt worden ist, wird im Schritt 106 gepriift, ob ein Abbruchkriterium erfiillt ist Dieses Ab- 
bruchkriterium kann z.B. eine vorbestiromte Anzahl an Wiederholungen der Schritte 103 bis 
10 105 sein oder eine kaum noch voihandene Anderung der Lage der Netzpunkte. Ist dieses 
Abbracbkriterium erfiillt, so endet das Segmentierungsverfahien im Schritt 107. Ist das 
Abbruchkriterium dagegen noch nicht erfiillt, so wird das Verfahien mit dem Schritt 103 
fortgefiihrt In anderen Ausfuhrungsfonnen konnte auch mit dem Schritt 104 fortgefahren 
werden. 

15 

Das nach dem Segmentierungsverfahren erhaltene Modell kann gespeichert und auf dem 
Monitor dargestellt werden. Es stellt die Segmentierung der in dem dreidimensionalen Daten- 
satz enthaltenen zu segmentierenden Struktur dar. 

20 In anderen Ausfuhrungsfonnen konnten die Schritte 104 und 105 mit dem Schritt 103 ver- 
tauscht werden. Dann wiirden zunachst die Kandidatenpunkte ermittelt werden mittels denen 
dann imterMinimierung dear extemen und der intemen Energie das Modell deforrrriert wird Im 
nachsten Schritt hatte ein Anwender die Moglichkeit, einzelne Dreiecksflachen bzw. einzdne 
Netzpunkte des deformierbaren Modells zu verschiebert Wenn die nachfolgende ttberpriifang 

25 des Abbnichkriteriums eatgibt, dass dieses Kriterium noch nicht erfullt ist, so wild nach Errnitt- 
lung neuer Kandidatenpunkte das Modell unter Minimising der Gesamteneargie, einschliefilich 
£nser» erneut def ormiert 
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Des weiteren kann in anderen Ausfuhrungsfonnen ein zusatzlicher Defoimationsschritt 105 
zwischen den Schritten 103 und 104 eingefiigt werden. Dies ermoglicht jeweils eine Kontrolle 
des Einflusses der Kandidatenpunkte und des Einflusses der verschobenen Dreiecksflachen 
auf die Deformation des Modells. 

5 

Die Vorziige des erfindungsgemaBen Verfahrens lassen sich anhand der Figuren 4 bis 7 er- 
klaien. In diesen Rguien ist eine zweidimensionale sagittale Schicht eines dreidimensionalen 
Datensatzes zu sehen, wobei in diesem Ausffihrungsbeispiel der Datensatz mittels eines Com- 
putertomographen erzeugt worden ist Ein defonnierbaies Modell M ist urn einen Riicken- 

10 wirbel 6 einer Wirbelsaule 7 angeoidnet Fig. 4 zeigt diese Anordnung vor Ausfuhrung eines 
Anpassungsverfahrens. Ein typisches Resultat der Deformation des Modells M mittels eines 
der in der Einleitung genannten Verfahien ist in Kg. 5 zu sehen. Im Bereich des iechten Ab- 
schnitts 8 des Modells M ist ein Segmentierungsfehler aufgetneten, da der Modellabschnitt 8 
auch Teile des benachbarten, unteien Ruckenwirbels 9 aufweist In Fig. 6 ist manuell versucht 

15 worden, den Segmentierungsfehler grob zu korrigieren, indem der Modellabschnitt 8 verformt 
bzw. zum Teil nach rechts verschoben worden ist Nach dieser groben Korrektur stimmt der 
Randbeieich 1 1 des Modellabschnitts 8 gut mit der Grenzflache der zu segmentierenden 
Struktur uberein, wohingegen der rechte Randbeieich. 13 des Modellabschnitts 8 weniger gut 
die Grenzflache annahert Fig. 7 zeigt das Modell M nach Durchfuhrung des etfindungsge- 

20 maBen Verfehrens, wobei die Gesamtenergie gemaB Gleichung (4) minimiert worden ist Der 
linke Randbeieich 11 ist kaum verandert woiden. Der rechte Randabschnitt 13 wurde da- 
gegen stalker deformiert und ist gut an die Grenzflache angepasst Die Ursadie fur die unter- 
schiedKche Behandlung der manuell veranderten Randabschnitte 11, 13 sind die Gewichtungs- 
faktoiCTt w t in Gleichung (8) bzw. (9). Der Randabschnitt 11 lag auf bzw. nahe der Grenz- 

25 flHche dear zu segmentierenden Struktur, so dass die Kantencfetektionsfunklion F (x, ) groB 
ist, was zu einan gtx>fien Gewichtungsfaktor fuhrt Der Randabschnitt 13 war dagegen weit 



-13- 



PHDE020307EP-P 



von der Grcnzflache entfemt, so dass F(x £ ) einen Meinen Wert liefert, was ednen Meinen 
Gewichtungsfaktor oder w t = 0 zur Folge hat Die Anpassung des Randbereichs 13 an die zu 
segmentierende Struktur erfolgte also allein durch Minirnierung der Energie + aJE m t. 

Dutch die Bewertung der manuell verschobenen Netzpunkte ist es nicht notwendig, jeden 
Netzpunkt bzw. jede Dreiecksflache des Modells moglichst genau auf die Oberflache der zu 
segmentierenden Struktur zu verschieben. Daher wird durch diese Bewertung eine Benutzer- 
korrektur mit geringem Zeitaufwand ennoglicht 
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BEZUGSZEICHENUSTE 

M deforniierbaresMcxJell 

1 BndvoHdbatun^- und Steuaprozessca: 

5 2 Spacher 

3 Glasfaffrkabd 

4 Monitor 

5 Tastatur 

6 Ruckmwrbel (zu segpnentienende Struktur) 
10 7 WirbeJsauie 

8 Abschnitt des defanmobaren Modells 

9 Rucken\\aifoeL 
11,13 RanddesAbschnitts8 
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PATENTANSPRUCHE 



1. VerfahiBi zum SegicnentieEin. aner dreidiirimaoaialen Struktur aus eirtem 

drediirienaorialQi, insbesondoB medizmischai Datensatz unter Boik^chtigMng^on 
Beniitzedcooekharm nit den Schritten: 

5 a) Berritstellen eines dreddimensionalen, deformierbaien Modells (M), dessen 
OberfLache durch ein Netz mit Maschen gebildet ist, die Netzpunkte an der 
Oberflache des Modells miteinander verbinden, 

b) Anordnen des Modells (M) an einer Stelle im dreidimensionalen Datensatz, an der 
10 sich. die zu segmentieiende Struktur (6) befindet, 

c) manuelles Verschieben von Netzpunkten, 

d) Neuberechnen der Netzpunkte des Modells (M) unter gewichteter Beriicksichtigung 
15 der manuell verschobenen Netzpunkte, 



20 
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Verfahimnach^^ 

- Bestimmenje dries 

beg^enzte Telflache, wobd jeder Kandidatenpunkt auf drier TdMachomomialen 
liegt, 

Zuordnen je dnes Gevvichtungpfaktois zu jedem verschcbenen Netepunkt, w^bdder 
Gevvichtungpfaktor, urn so gpofier ist, je Werner der island des voschobmen 
Netqxnikteszu drier Gre^^ 

NeubeiadinmderNeb^^ 

Kandidatenpunkte, der verschobenen Nefepunkte und dor zugsordneten 
GevAichhin^aktoiea 



Verfaihtenrradh dadiKch gekenii2d^ daasimSchritt d) dieNet2punkte 

unter Miriiinierung einer gewLchtetai Summe aus externa: imdmtemer Energie und 
duo: Energfe die die manuell Netzpunkte berucikachtigt, neubeoechnet 
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4. Bildvetaibatungg^^ 

- einer Speichereinheit zum Speichem eines deformierbaien Modells, dessert 
5 Oberflache durch ein Netz mitMaschen gebildet ist, die Netepunkte an der 

Oberflache des Modells miteinander verbinden, und zum Speichem eines 
dreidimensionalen Datensatzes, insbesondere eines medizMs chen Datensatzes, 

- einer BilcMedo^ibeanheit zur Wiedergabe eLner zu segrcntieDendenbzw 
9egpientierten Struktur und des defomuabaren ModeJls> 

10 - ana: Bo^chniin^anhat zur Neuberednnung der Netzpunkte des Modells unter 
ggwichteter Boikksidttigpng \un manudl veschobenm Nefepunlderv 

- einarArordnung^e^ 

c^dimenaonalen Datensafc^ an der ach die zu segpnenliecende Struktur b^ndet, 

- oner Steuerenhat zur Steuerung der Spacteenhat, da: Bildvviedergpbeeinhat, der 
15 Ba^chmin^arihat und da: Anoixkixin^anhat antqprechend dai f olgpnden 

Schritten; 

a) BereitstdLen eines c^dimenaanalai, defbmuerbaren Modells dessen Oberflache 
dutch ein Netz mit Maschen g&ildet ist, die Net2punkte an der Oberflache des 

20 Modells imteinanda^ 

b) Au»x3nai des Moddls an einer St^ 
die zu segpnentierende Struktur befindet, 

c) manudles Verschidsen \on Netzpunkten, 

d) NeubetechnenderNetepu^ 

25 manudl TCrschobenen Netzpunkte 
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Computerpiograrnrn fur eine Steuereinheit zur Steuerung einer Speichereinheit, einer 
Bildwiedergabeeinheit, einer Berechnungseinheit und einer Anordnungseinheit einer 
BMveambdttmgseinriditung zur Duichfuhrung des Verfahrens nach Anspruch 1 
entsprechend den folgenden Schritten: 

a) BCTatstellen dues cteidiniaiskHTaler^ (Mormiefcarea MorielK dessen Oberflache 
dutch ein Netzmit Maschen gEbildet ist, die Netzpunkte an dec Oberflarhe des 
Modells iruteinander \etrindetv 

b) Anordnen des Modells an drier Stdle im. drdcKmensioraLen D aiensatz, ander sich 
die zu segpcieniierende Strukturbeflndet, 

c) manuelies Veischieben \cnn Netzpunkten, 

d) Neuberechnea der Nefczpunkte des Modells unter g2v\ichteter Berucicsichtigung der 
manueJl voschobeno^ Nefepunkten 
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Vecfehren zum Segpnentiecen ener drddimensianalen Struktur aus etnem dradimenaonaleri 
Datensatz unter BerudssicMgungvon Benutzakoriektuim 

5 Die Erfindung betrifft ein Verfehren zum Segmentieren einer dreidimensionalen Struktur aus 
einem dreidimensionalen, insbesondere medizinischen Datensatz unter Beriickachtigung von 
Benutzeikorrekturen. Das Verfahren wird mit Hilfe eines dreidimensionalen deformierbaren 
Modells, dessen Obetflache durch ein Netz mit Netzpunkten und diese verbindende Maschen 
gebildet ist, durchgefuhrt Nachdem das Modell an der Stelle im dreidimensionalen Datensatz, 

10 an der sich die zu segmentierende Struktur befindet, angeordnet und ggf. Positionen von 
Netzpunkten durch bekannte Segmentieningsverfahren geandert worden sind, konnen 
beliebige Netzpunkte manuell verschoben werden. Die Netzpunkte des Modells werden unter 
gewichteter BeriicksicMgung der manuell verschobenen Netzpunkte neu berechnet 
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Fig. 4 
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